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Abstract. We describe an experimental apparatus capable of achieving a high loading rate of strontium 
atoms in a magneto-optical trap operating in a high vacuum environment. A key innovation of this setup 
is a two dimensional magneto-optical trap deflector located after a Zeeman slower. We find a loading rate 
of 6 X 10® s“^ whereas the lifetime of the magnetically trapped atoms in the ®P 2 state is 54 s. 

PACS. 3 7.10.De, 37.10.Gh 


1 Introduction 

Cooling and trapping of alkaline-earth atoms, and in the 
particular case of strontium atoms, have been rapidly de¬ 
veloped over the last two decades [T1IU[31I11[51[S1[71[H1IS] • These 
two valence electrons atoms have a ground state with zero 
electronic spin and weakly allowed singlet-triplet transi¬ 
tions. They offer interesting alternatives to the more com¬ 
monly used alkali atoms and open new fields of research 
for cold and ultracold gases. In the field of ultra high pre¬ 
cision spectroscopy, tremendous progress has been made 
recently to reduce the frequency instabilities and system¬ 
atic shifts of atomic clocks. Currently, the best results are 
obtained with strontium lattice clocks (see nnun] for first 
implementations and for reviews) where the ac¬ 

curacy and clock frequency difference were measured at 
the 10“^® level [TKlllfij . Those developments are of im¬ 
portance not only for time keeping, but also for preci¬ 
sion tests in fundamental physics HZlllHllll]. Using dif¬ 
ferent strontium isotopes, test of the equivalence princi¬ 
ple and search for spin-gravity coupling effects have been 
reported |2nj- In studies related to quantum gases, pro¬ 
duction of a ®^Sr Bose-Einstein condensate (BEC) was 
achieved by the Innsbruck and the Rice groups [5|[5]. These 
were rapidly followed by the production of quantum gases 
of different Sr isotopes and their mixtures pni^lMl[24] . 
Bosonic strontium atoms are textbook examples of scalar 
BEC, whereas the ®^Sr Fermion has a large 7 = 9/2 nu¬ 
clear spin that allows the study of Kondo physics [25] . 
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SU(N) magnetism [5^ and gauge field [?7ll^ . The lack 
of electronic and nuclear spin in the ground state of the 
bosonic isotopes reduces the resonant photon scattering 
phenomena to its simplest component: Rayleigh scatter¬ 
ing (ie., no Raman scattering). In this context, it has been 
shown that the photon coherence is preserved during mul¬ 
tiple photon scattering [29]. Furthermore, it is destroyed 
only by the intensity saturation effect (SO]- Finally, the 
narrow intercombination line of ®®Sr was used to reveal 
cooperativity and (super)flash effects in forward scatter¬ 
ing of an optically thick cold cloud mm- 

In this article, we describe a new experimental ap¬ 
paratus that achieves a high loading rate of a strontium 
magneto-optical trap (MOT) in a high vacuum environ¬ 
ment. Indeed, we measure a loading rate of 6 x 10® s“^ 
which is, to the best of our knowledge, the highest rate 
reported for a strontium setup. Such a high rate consider¬ 
ably reduces the loading time of the MOT. This is helpful 
to decrease the duty cycle in optical lattice clocks, mini¬ 
mizing the Dick effect (see for example [33] and references 
therein), or to produce degenerate quantum gases with 
higher number of atoms |53|. High loading rates also lead 
to samples that are optically thicker, which are important 
to explore cooperative phenomena [32] . 


We describe the different components of our experi¬ 
mental apparatus in section [^ One key feature is a two 
dimensional MOT deflector located after a Zeeman slower 


(sections |2.3| and 2.4). It increases the brightness of the 
slow beam of atoms and deflects it by an angle of 30°. The 
beam is redirected into a differential pum ping tube before 
it enters into the science chamber (section 2.51. This setup 
allows us to maintain a low residual background pressure 
at the position of the MOT, resulting in a 54 s lifetime of 
the ®P 2 state magnetic quadrupole trap. The details of the 
laser systems addressing the various transitions (shown in 
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Fig. 0 are given in section In particular, the laser sys¬ 
tem, operating on the —> ^Pi dipole allowed transi¬ 

tion at 461 nm, is locked on a separate Sr atomic beam. 
We tune the frequency lock point by varying the magnetic 
field. This allows us to address the different strontium iso¬ 


topes (section 3.1). Finally, we modulate the phase of both 
the 707 nm and 679 nm repumping lasers. We generate fre¬ 
quency sidebands which address, at the same time, ®®Sr 
and the complete hyperfine structure of ®^Sr. This aspect 
is discussed in section 13.31 


2D OM before entering the Zeeman slower, where its longi¬ 
tudinal velocity is reduced to 60 m/s. This value is within 
the velocity capture range of the MOT. Finally, the atoms 
are deflected and guided into the MOT chamber by a 2D 
MOT. All these preparatory stages, described in detail in 
the following sections, use the broad blue dipole allowed 
transition. 
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Fig. 1. Energy levels and transitions of interest of strontium. 
The excited states decay time and wavelength are indicated on 
the graph. 


Fig. 2. Three dimensional schematic view of the vacuum ap¬ 
paratus. 


Excluding the strontium oven, the other elements of 
our apparatus are analyzed and optimized in terms of the 
atom loading rate Lq into the MOT. Lq is extracted from 
measurements of the fluorescence signal of the cold gas. 
For this purpose, we use a resonant 10 mW blue probe 
beam with a waist of 3.5 cm that is larger than the cold 
cloud size. From the fluorescence signal, we deduce the 
number of cold atoms in the MOT. To avoid systematic 
errors in the atom counting, we calibrate the fluorescence 
signal with the absorption imaging, by a CCD camera, of a 
cold cloud containing a small number of atoms, using the 
same probe. The loading rate values are extracted from 
a fluorescence measurement, done 5 ms after turning on 
the MOT beams. Since this time is much shorter than the 


2 Experimental apparatus 

The relevant energy levels and transitions for laser cool¬ 
ing and trapping of strontium atoms are shown in Fig. 

From the singlet ground state ^Sg, there are two cool¬ 
ing transitions. The first one is the blue dipole allowed 
transition ^Sg —?► ^Pi at Ag = 461 nm. Its linewidth is 

Fb — 27r X 32 MHz and the saturation intensify is 42.5 mWcnjij^g^j^ Under these conditions, we do not expect 


characteristic MOT lifetime (section 2.5), the number of 
atoms in the MOT increases linearly. Moreover, the max¬ 
imum number of atoms in the MOT during the loading 
ra^ measurement is 3 x 10^ with a density that is lower 


This transition is u s ed t o m anipu late the effusive atomic 
beam (sections 2.2 2.3 and 2.4). It is also used in the 
first stag e of laser cooling in a magneto-optical trap (sec¬ 
tion [2^. The other cooling transition is the red intercom¬ 
bination line ^Sg —>■ ^Pi at Ar = 689 nm. Its linewidth 
Fr = 27r X 7.5 kHz is much narrower than the blue tran¬ 
sition. Its saturation intensity 3 pWcm”^ is also much 
weaker. This transition is used in the final cooling stage 
for both the bosonic ®®Sr and the fermionic ®^Sr isotopes 
(section |2.6[ ). 

A sketch of the experimental setup is shown in Fig. 

It consists of an oven and four stages of cooling and ma¬ 
nipulation of the atomic beam, namely, a two dimensional 
optical molasses (2D OM), a Zeeman slower, a two di¬ 
mensional magneto-optical trap (2D MOT), and finally a 
magneto-optical trap (MOT). The effusive atomic beam, 
extracted from the oven, is first transversely cooled by the 


any bias on the atom counting due to multiple photon 
scattering [55] . 


2.1 Strontium oven 

The oven is made of a stainless-steel cylinder heated at 
a maximal temperature of 550 °C. This temperature is 
stabilized to ensure a constant flux of the outgoing effu¬ 
sive atomic beam. The beam is generated in a multichan¬ 
nel nozzle made of nine hundred Monel400 micro tubes 
clamped in a 1 x 1 cm square area. The length of each 
tube is 1 cm with inner and outer diameters of 200 /im 
and 300 fim respectively. We keep the nozzle at a tem¬ 
perature of roughly 50 °C above the oven temperature to 
avoid clogging the tubes. 





















T. Yang et al.: A high flux source of cold strontium atoms 


3 


In Fig. [ga), we give two examples of the transverse 
transmission spectra of the atomic beam. These spectra 
are taken at 10 cm after the oven nozzle. They reflect the 
transverse velocity distribution of the atomic beam. The 
velocity distributions have a triangular-like shape, a char¬ 
acteristic of the collisionless regime inside the nozzle micro 
tubes [30] • As a result, we observe that the brightness of 
the atomic beam increases exponentially in the tempera¬ 
ture range of interest [Fig. [^b)]. 




Fig. 3. (a) Transverse spectra of the atomic beams after the 
oven. The solid red (blue) curve corresponds to a tempera¬ 
ture of T = 490 °C (T = 550 °C). To compare both profiles, 
we rescale the low temperature curve such that both spectra 
show the same amplitude. This rescaling operation is repre¬ 
sented by the dashed red curve. The two vertical solid lines 
indicate the position of the ®®Sr and ®®Sr resonances at 0 MHz 
and —125 MHz respectively, (b) Atomic beam brightness as a 
function of the oven temperature. 


2.2 Two dimensional optical molasses 

The effusive atomic beam, produced in the oven as ex¬ 
plained in the previous section, is transversely cooled by a 
two dimensional optical molasses (2D OM) located 15 cm 
after the nozzle. We use two pairs of retroreflected beams 
having the same elliptical shape (beam waists: 1 cm and 
0.5 cm). The major axis of the ellipse is oriented along 
the atomic beam. The performance of the 2D OM cool¬ 
ing is measured in terms of the atom loading rate Lg into 
the MOT. This quantity is plotted in Fig. for differ¬ 
ent values of the frequency detuning and the total power 
of the cooling beams. The optimum performance is found 
for a detuning of —18 MHz = —0.55/6. It corresponds to a 
three-fold improvement of the MOT loading rate. We note 
that further improvements may be expected at a higher 
laser power. Here, we are limited by the total available 
laser power (section [3d^. 


After the 2D OM, the transversely cooled atoms pass 
through a 10 cm long tube with an inner diameter of 1 cm. 
The main purpose of this tube is to block the unwanted 
atoms with high transverse velocities. By this means, we 
also create a differential pressure environment between the 
oven and the subsequent stages of the vacuum apparatus. 



Fig. 4. MOT loading rate as a function of the 2D OM total 
optical power for different detunings indicated on the graph. 
Since the beams are retro-reflected, the total power is twice 
the incident power. 


2.3 Zeeman slower 

The longitudinal velocity of the atomic beam is reduced 
using a 32 cm long Zeeman slower (see Fig.[^. The inho¬ 
mogeneous longitudinal magnetic field profile of the Zee- 
man slower is created by winding flat copper wires in a 
tapered way. At the optimum electrical current of 6.6 A, 
the magnetic held of the Zeeman slower varies from 500 G 
at the entrance to —100 G at the output. The optical beam 
of the Zeeman slower has a frequency detuning of —13/1, 
and a power of ~ 110 mW. The laser beam enters into the 
apparatus through an optical grade z-cut sapphire view¬ 
port. We permanently heat the viewport at a temperature 
of 375 °C to prevent the deposition of Sr atoms. 

We scan the frequency of a probe around the blue res¬ 
onance to measure the longitudinal velocity distribution 
of the atomic beam at the output of the Zeeman slower. 
The probe propagation axis is at an angle of 70° with re¬ 
spect to the atomic beam, limiting the velocity resolution 
to ~ idr’b/Zcb. Here, T'b/fcb = 15 m/s, and kb = 27r/A6 
is the wavevector of the probe. Some examples of the ve¬ 
locity distribution are given in Fig. [^a). When the Zee- 
man slower is turned off, the velocity distribution is in 
good agreement with the expected longitudinal distribu¬ 
tion calculated for an oven at T = 450 °C. Once the Zee- 
man slower is turned on, the velocity of a large fraction 
of the atoms is brought down, contributing to a peak at 
^ IFb/kb [see the blue curve in Fig. [^a)]. We note that 
the Zeeman slower is effective for vGocities lower than 
~ db/ft/fcb. This value is in agreement with the Zeeman 
shift imposed by the magnetic field at the entrance. 

In Fig.[5];b), we show in detail the performance of the 
Zeeman slower by plotting the loading rate of the MOT 
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(a) 




Fig. 5. (a) Atomic longitudinal velocity distribution measured 
after the Zeeman slower. We use a dimensionless velocity unit 
where Fb/kb = 15 m/s. The red solid line is the atomic veloc¬ 
ity distribution without the Zeeman slower optical beam. The 
black dashed line is the theoretical distribution for an oven at 
T = 450°C. The blue line is the velocity distribution with the 
Zeeman slower operating optimally, (b) MOT loading rate as 
function of the Zeeman laser power when the magnetic held is 
turned on (red dots and solid lines) and when it is turned off 
(blue stars and solid lines). 

as a function of the optical power of the Zeeman slower. 
Two situations are analyzed; when the Zeeman slower 
magnetic field is off (blue stars) and when an optimum 
magnetic field prohle is applied (red dots). At low power, 
i.e. < 30 mW, the two situations are equivalent, meaning 
that the weak radiation pressure force is not enough to 
impose a deceleration of the atoms compatible with the 
variation of the magnetic field inside the Zeeman slower. 
The Zeeman slower is thus inoperative. Above a power of 
60 mW, the Zeeman slower reaches its optimum perfor¬ 
mance, showing a 20-fold improvement of the MOT load¬ 
ing rate compared to the situation where the magnetic 
field is turned off. Overall, the Zeeman slower improves 
the loading rate in the MOT by a factor of 300. 

The magnetic field of the Zeeman slower is confined 
at its vicinity using f-metal sheets (see Fig.[^. This mag¬ 
netic field shielding box does not affect the performance of 
the Zeeman slower. Rather, it prevents a distortion of the 
magnetic field lines of the two dimensional atomic beam 
deflector located just after the Zeeman slower. 

2.4 Two dimensional magneto-optical trap and beam 
deflector 

After the Zeeman slower, the atoms enter a two dimen¬ 
sional magneto-optical trap (2D MOT). This device aims 


to transversely cool, focus and guide the atomic beam 
along the zero magnetic field line located along its sym¬ 
metry axis To deflect the atomic beam, the 2D 

MOT axis is purposely kept at an angle of 30° off the Zee- 
man slower axis , in the horizontal plane (see Fig. [^and 

t . At the output of the 2D MOT, the atomic beam passes 
rough a differential pumping tube with an entrance ra¬ 
dius of 2.5 mm and an acceptance angle of 10 mrad. 

Under this original arrangement, our 2D MOT has 
at least four advantages. First, by deflecting the atomic 
beam, we prevent the intense Zeeman optical beam to 
cross, and thus to disturb, the MOT in the science cham¬ 
ber. Second, the uncooled atoms in the atomic beam, in¬ 
cluding the unwanted Sr isotope, do not enter the science 
chamber. Combined with the differential pumping tube, it 
allows us to isolate the principal part of our experimental 
setup, namely the science chamber, from the rest of the 
vacuum apparatus. Hence, we get an ultra high vacuum 
environment in the science chamber. This allows us to 
achieve a long lifetime of our cold sample for future evap¬ 
orative cooling towards degenerate gases. Third, we gain 
an additional optical axis in the science chamber along 
the atomic beam. Finally, once the MOT is loaded, the 
2D MOT beams can be turned off so that the atomic flux, 
entering into the science chamber, drops to zero. Hence, 
this design ensures a high atomic flux, but does not require 
any extra mechanical blocker for the atomic beam. 



Fig. 6. Schematic view of the 2D MOT atomic beam deflector. 
The dashed line gives the direction of the initial atomic veloc¬ 
ity. Two pairs of coils generate a 2D quadrupole field with a 
line of zero magnetic field along the 2 :-axis a,t a. 6 = 30° angle 
with the incident atomic beam. The z-axis is also the outgoing 
direction from the 2D MOT. The two pairs of the retroflected 
MOT beams - a~ polarizations) are perpendicular to each 
other in the x — y plane. After the 2D MOT, the deflected 
atomic beam passes through a differential pumping tube. 

A schematic view of the 2D MOT is given in Fig. 
The two dimensional quadrupole magnetic held is gener¬ 
ated using four square shaped coils. The inner and outer 
sizes of the coil are 15 cm and 19 cm respectively. Each coil 
is formed by 22 turns x 20 layers of 1 mm diameter cop- 








T. Yang et al.: A high flux source of cold strontium atoms 


5 


per wire. The distance between each parallel pair of coils 
is 20 cm. This configuration gives a homogenous magnetic 
field gradient of dB/dx = 5 G/cm in the interacting zone 
for 1 A of current. The cooling beams have an elliptical 
shape with an aspect ratio of 3:1. The major axis is ori¬ 
ented along the atomic beam with a waist of 3 cm. Using 
retro-reflected beams, the total available optical power is 
220 mW. After optimization, we found that 60% of the 
total power has to be fed into the horizontal beam pair. 
Indeed, along this axis, larger radiation pressure forces are 
required to deflect the beam. The remaining power is sent 
into the vertical beam pair, where only the cooling and 
focusing of the atomic beam are required. 

As we did for the 2D OM and the Zeeman slower, we 
analyze the performances of our 2D MOT beam deflector 
in terms of Lq, the loading rate in the MOT. In Fig.j^a), 
we show Lq as a function of the total 2D MOT optical 
power for different detuning. The magnetic held gradient 
is set to the optimal value of 7.5 G/cm. As a general fea¬ 
ture, we observe a saturation of the loading rate above 
an optical power of ~ 140 mW. Surprisingly, it occurs 
at a modest saturation intensity parameter of 0.16. This 
value is calculated for a frequency detuning of —0.971,. 
We attribute this saturation behavior of the loading rate 
to a loss mechanism induced by spontaneous decay from 
the ^Pi state to the underlying long-lived state. This 
optical pumping mechanism is known to reduce the blue 
MOT lifetime. In principle, this can be circumvented using 
repumping lasers, as it is successful ly d one for the MOT 
in the science chamber (see section 2.5 for more details). 


Unfortunately, it is inefficient here since the dwell time of 
the atoms in the 2D MOT is shorter than the time for 
repumping. 


We perform numerical simulations to confirm the im¬ 
portant role of the ^Pi —> ^D 2 spontaneous decay as a loss 
mechanism in the 2D MOT. Starting from realistic veloc¬ 
ity and spatial probability distributions at the input of the 
2D MOT, we compute the classical trajectory of an atom 
subjected to the following mean radiation pressure force: 
F = {F^j^ + F^_)x -I- (7^y+ + Fy_)y. The component of 
the force along Ox axis is 


^ _ ^ nkFi, _ I^{y,z)/Is _ 

2 I + {2{St kvxT k‘b'x/h)/rb)^’ 


( 1 ) 


where &' is the magnetic field gradient and Ix{y,z) the 
intensity of the corresponding beam. One has a similar 
expression for the components of the force along the Oy 
axis, Fy±. We stop the simulation once the atom reaches 
entrance plane of the differential pumping tube. A suc¬ 
cessful transmission count fulfills two criteria. First, the 
atom must be loc ated at t he tube entrance within an ac¬ 
ceptance angle, + Vy/v^, smaller than 10 mrad. Sec¬ 
ond, the atom should not have spontaneously decayed to 
the ^D 2 state. For this purpose we evaluate the probabil¬ 
ity of a decay using the transition rate npFo, where tfp 
is the averaged population of the excited ^Pi state and 
Fd = 4 X 10 ^ s“^ is the ^Pi —D 2 decay rate. Gonsider- 
ing the transmission of the atoms through the 2D MOT 


with decay to the ^D 2 state included, the numerical sim¬ 
ulation is in qualitative agreement with the experiment; 
it shows a saturation plateau for an optical power above 
^ 140 mW [see the black open circles and the solid lines 
in Fig. [^b)]. If the decay to the ^D 2 state is removed 
from the simulation, we observe a qualitatively different 
behavior, namely a larger transmission with a further in¬ 
crease at high optical power [see the blue open stars and 
the dashed lines in Fig. [^b)]. 


(a) 




Fig. 7. (a) Loading rate of the MOT as a function of the total 
optical power of the 2D MOT for different frequency detuning 
indicated on the graph. Since the beams are retro-reflected, the 
total power is twice the incident power. The 2D MOT mag¬ 
netic field gradient is 7.5 G/cm. (b) Numerical simulation of 
the 2D MOT transmission percentage for a frequency detuning 
of —0.971,. The total number of atoms used in the simulation is 
10®. The black open circles and the solid lines give the results 
considering optical pumping to the long lived state ^D 2 . The 
case without optical pumping corresponds to the blue open 
stars and the dashed lines. 


2.5 Blue magneto-optical trap 

The atomic beam is transversely cooled, focused and de¬ 
flected by the 2D MOT as described in section |2.4[ After 
this stage, the atomic beam goes into the science cham¬ 
ber with a typical longitudinal mean velocity of 471,/fcb ~ 
60 m/s. This velocity is compatible with the velocity cap¬ 
ture range of our MOT operating on the blue dipole al¬ 
lowed transition. It is, however, too high to be directly 
captured by the red MOT operating on the intercombina¬ 
tion line (the typical capture velocity of the red MOT is 
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1 m/s). Thus, the atoms are initially cooled and trapped in 
the blue MOT before being transferred into the red MOT 
for further gain in the phase-space density. The operation 
of the red MOT is described in the section 12.61 


Our blue MOT setup is made of three pairs of retroflected 
beams with a waist of 1.1 cm, and a maximum total inten¬ 
sity of 76 mWcm“^. The typical performances of the blue 
MOT, in terms of the loading rate, are given in Fig.|^ We 
find a maximum atom loading rate of about 6 x Im s“^ 
at a magnetic field gradient of 37.5 G/cm and a frequency 
detuning of —58 MHz. 


The lifetime of the cold atomic cloud depends strongly 
on the presence or absence of repumping lasers. Indeed, 
the ^Sq —Pi cooling transition is not perfectly closed 
because of the presence of an underlying ^D 2 state. As 
a result, atoms are lost by spontaneous decay into the 
long-lived triplet state ^P 2 , with a typical MOT loss rate 
of 50 s“^. This value strongly depends on the ^Pi state 
population, and therefore on the MOT beams intensity 
and frequency detuning. With repum ping lasers (a de¬ 
tailed description is given in the section 3.3), this loss rate 
is reduced. An example of the temporal evolution of the 
atoms number in the MOT, after turning on the MOT 
beams, is shown in Fig. [^a). Fitting the data with an 
exponential function, we und a loss rate of 14 s“^, lim¬ 
iting the atoms number in the MOT, at the stationary 
regime, to iV ~ 5 x 10®. Such a high value of the loss 
rate is rather surprising. We believe it could be explained 
by unwanted escape channels in the MOT or incomplete 
repumping of the atoms from the metastable states. We 
shall exclude losses due to light assisted collisions since the 
two-body collisional loss coefficient reported in the litera¬ 
ture is /3 ~ 2 X 10“^° cm® /s gnus]. With a peak spatial 
density of n ~ 3 x 10® cm“®, this leads to a small light as¬ 
sisted collision loss rate, estimated to be /3n~ 0.6 s“®. We 
exclude the background collision losses as well. This is in¬ 
ferred from the measurement of the atomic cloud lifetime 
in the ®P 2 state trapped in the magnetic quadrupole field. 
This experiment is done by first loading the MOT without 
repumping lasers for 1 s. After this loading sequence, we 
hold the atoms for a variable time before turning on the 
repumping laser. The population in the magnetic trap is 
then transferred to the ground state and measured in the 
MOT. We find an atoms loss rate of 18 x 10“® s“^ in the 
magnetic trap [Fig. [^b)]. This loss rate is in agreement 
with loss mechanisms induced by the blackbody radiation 
at room temperature through the ®P 2 —> ®Di, 2,3 transi¬ 
tions g2]. This gives an upper limit to our background 
collision loss rate and confirms the good vacuum environ¬ 
ment in the science chamber. 


Finally, we note that higher atoms number can be ob¬ 
tained using an atom shelving technique by switching off 
the repumping lasers. In this technique, the atoms are 
loaded in the long lived ®P 2 magnetic trap rather than in 
the blue MOT. g^l^l44] . 


(a) 




Fig. 8. The atomic loading rate in the MOT as a function of 
MOT beam total power. Since the beams are retro-reflected, 
the total power is twice the incident power, (a) Each curve 
corresponds to a different frequency detuning indicated on the 
graph. The magnetic field gradient is 37.5 G/cm (b) Each curve 
corresponds to a different MOT magnetic field gradient indi¬ 
cated on the graph. The frequency detuning of the MOT beams 
is — l.STf,. 


2.6 Red magneto-optical trap 


The final laser cooling sequence is performed in a red 
MOT operating on the ®So —>■ ®Pi intercombination line 
at 689 nm. We use similar techniques to those reported in 
the early studies by the JILA and the Tokyo groups [TJ 
H5] . In preparation for the transfer into the red MOT, the 
temperature of the cold gas in the blue MOT is reduced by 
decreasing the power of the blue MOT cooling beam down 
to 5 mW in 50 ms. The laser power is then kept at this 
value for 100 ms (see the time sequence in Fig. 10). The 
final temperature in the blue MOT is about 2 mK/corre¬ 
sponding to an rms velocity of D « 0.5 m/s |46] . In the red 
MOT, the minimal stopping distance is u^/(2|a|) « 1 mm, 
where jaj = hkr/{2rr) « 155 ms“® is the acceleration in¬ 
duced by the maximal radiation pressure force. This sim¬ 
ple calculation shows that, first, the stopping distance can 
be shorter than the typical laser beam size (1.4 cm in our 
case). Second, it is crucial to maintain a large radiation 
pressure force across the initial velocity distribution of the 
cloud. For this purpose, the laser frequency linewidth is ar¬ 
tificially broadened to 2 MHz« 3kv by frequency modula¬ 
tion m- The total power of the red MOT beams is 30 mW 
and the modulation rate is 25 kHz. Thus, the intensity 
per red MOT beam and on each comb is 120 p,Wcm“®, 
which is 40 times the saturation intensity. This configu¬ 
ration ensures a large radiation pressure force during the 
transfer into the red MOT. Finally, the artificial broad¬ 
ening is turned off to optimize the temperature and/or 
the spatial density. In the following section, we discuss 
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(a) 




Fig. 9. (a) Temporal evolution of the number of ®®Sr atoms 
in the blue MOT. The MOT beams are turned on at time t = 
0. (b) Temporal evolution of the atoms number magnetically 
trapped in the ®P 2 state after 1 s of loading. The 2D MOT 
beams are turned off at t = —100 ms. For both graphs, the 
dots are the experimental values and the solid red lines are 
exponential fits of the data. For both graphs, the magnetic 
field gradient is set at 45 G/cm. 


in further details, the operation of the red MOT for the 
®®Sr and ®^Sr isotopes. The difference between these two 
cases comes from the presence of a Zeeman state man¬ 
ifold and a hyperfine structure in the excited state for 
the fermionic isotope, introducing more complexity in the 
cooling scheme. 
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Fig. 10. Time sequence for the loading and cooling of ®®Sr in 
the red MOT. 


2.6.1 The ®®Sr boson 

The detailed time sequence of the red MOT for the ®®Sr 
bosonic case is shown in Fig. At the final stage of 
the broadband cooling, we typically transfer 50 percent 
of atoms into the red MOT. Then, the red MOT laser is 


switched to a single frequency with a detuning of —200 kHz 
« —660).. The MOT magnetic field gradient is ramped up 
from 0.75 G/cm to 4 G/cm in 100 ms to compress the cold 
cloud. At the final stage of the cooling sequence, we find 
3.5x10® atoms in a volume of 0.15 mm®, at the tempera¬ 
ture of r = 3 p-K. The spatial density is n = 2.5 x 10®^ cm“®, 
leading to a phase space density of 3 x 10“®. Finally, we 
find a density dependance of the MOT temperature with a 
slope of dr /An ~ 6 x 10“®^ cm®pK. This shows that extra 
heating due to collective effects, such as multiple photon 
scattering, are still in play, even if the narrow transition 
has a linewidth in the same order of magnitude as the 
recoil frequency w,., namely T). = O.Ow^. 


2.6.2 The ®'^Sr fermion 


In contrast to the bosonic isotope, the ®’^Sr fermion has a 
nuclear spin 7 = 9/2 which induces Zeeman sublevel mani¬ 
folds and hyperfine splitting of the excited state. As a con¬ 
sequence, the MOT is less efficient and technically more 
challenging. We use the transition F/ = 9/2 —>■ F/ = 11/2 
to provide the necessary magneto-optical trapping force. 
However, because of the narrow resonance linewidth to¬ 
gether with the presence of the MOT magnetic field, op¬ 
tical pumping among the Zeeman substates brings the 
atoms out of resonance. As initially proposed in Ref. [2], 
we ensure proper remixing of the Zeeman substate by 
adding a stirring laser tuned on the transition Fg = 9/2 —)■ 
Ff.=%l2. The stirring laser follows the same optical path 
as the MOT beams, thus it also participates in laser cool¬ 
ing. 


Using acousto-optic modulators (AOMs), we control 
independently the detuning and th e in tensity of the stir¬ 
ring and MOT lasers (see section 3.2 for more details). 
The loading and cooling sequence of the ®^Sr red MOT 
is sketched in Fig. After a blue MOT phase iden¬ 
tical to that for the bosonic isotope, we first decelerate 
the atoms in an optical molasses for 10 ms. Both stirring 
and MOT lasers are frequency broadened to 2 MHz and 
brought to their maximum power, i.e. 25 mW and 12 mW 
respectively. After the optical molasses stage, we apply 
a magnetic field gradient of 2 G/cm for 40 ms to com¬ 
press the atomic cloud. Then, we reduce the laser pow¬ 
ers to half of their initial values in 30 ms in order to de¬ 
crease the cloud temperature before switching to the single 
frequency operation. During this final cooling stage, the 
frequency detunings of the stirring and MOT lasers are 
—350 kHz and —220 kHz respectively, and their powers 
are progressively reduced to 3 mW and 2 mW respectively. 
We get 3x 10^ atoms at temperature of 6 pK in a volume 
of ^ 1 mm®. The reduction in the atoms number com¬ 
pared to the ®®Sr case is consistent with the difference in 
the natural isotopic abundance. However, the larger cloud 
volume and the higher temperature indicate that cooling 
and trapping are less efficient for the fermionic isotope. 






















T. Yang et al.: A high flux source of cold strontium atoms 


Blue MOT 


-2.2 MHz 
-0.2 MHz f 
Power, 

9/2 (raW) - 
Power, 
Il/2(mW)- 

B (G/cm) - 


Blue 
power (mW) 


75 


r 


-100 


Low power 
blue MOT 


45 


Broadband 
red MOT 


Single band 
red MOT 


25 


12 


0.75 


Final stage 
red MOT 


100 200 300 

Time (ms) 


400 


500 


Fig. 11. Time sequence for the loading and cooling of ®^Sr in 
the red MOT. 


3 Laser systems 


This section is devoted to the description of the laser sys¬ 
tems. The 461 nm laser, addressing the ^Sq —^Pi tran¬ 
sition, is composed of a master laser at 922 nm, a series 
of low and high power infrared (IR) optical amplihers and 
two second harmonic generation (SHG) resonant cavities. 
The blue master laser is frequency locked on a ®®Sr atomic 
beam. The lock point can be tuned via Zeeman effect using 
a pair of Helmholtz coils. This design allows us not only 
to change the global detuning of the laser, but also to ad¬ 
dress the different strontium isotopes (section 3.11. The 
high frequency stability required for the -T '^Pi in¬ 
tercombination line at 689 nm prevents us to implement a 
similar Zeeman based frequency tuning. Here, the different 
transitions of interest for ®®Sr and ®^Sr are addressed us¬ 
ing AOMs and a 1.2 GHz electro-optic modulator (EOM, 
see section 3.21. The red master laser is frequency locked 


on a high-hnesse cavity. Its long term drift is compen¬ 
sated using a saturated fluorescence spectroscopy setup 
on an atomic beam. Finally, two lasers at 707 nm and 
679 nm are tuned on the ^P 2 —>■ ^Si and ^Pg —>■ ^Si tran¬ 
sitions respectively for metastable state repumping during 
the blue MOT phase. With a single electro-optic modu¬ 
lator operating at high modulation index, we are able to 
address all the hyperfine sublevels of the ®^Sr ^P 2 state 
(section |3.2[ ). 


3.1 Blue laser 

The optical setup of the blue laser system is shown in 
Fig. Five output ports address the different parts of 
the experimental apparatus. AOMs and mechanical shut¬ 
ters are used to turn on and turn off the lasers beams. The 
primary laser is a commercial solid-state laser (Toptica 
TA-SHG Pro) which provides 400 mW output at 461 nm. 
It is composed of a 922 nm external cavity diode laser 
(FCDL) used as a master laser, a 1.5 W IR tapered am- 
pliher (TA), and a SHG resonant cavity. A 3 mW IR out¬ 
put port of the primary laser is used to inject a 30 mW 
slave laser, after a single pass through two AOMs respec¬ 
tively at —45 MHz and —90 MHz. This laser beam injects 


a 1 W TA after a 60 dB isolator. Then, the beam is fre¬ 
quency doubled using a home-made bow type ring cavity 
containing a 20 mm long periodically poled potassium ti- 
tanyl phosphate (PPKTP) crystal. We get 200 mW power 
at 461 nm. This beam, used for the Zeeman slower, is red 
detuned by —2 x 135 = —270 MHz with respect to the 
primary blue laser. Overall, the Zeeman slower beam is 
detuned by—408 MHz with respect to the frequency lock 
point. 


The master laser is frequency locked using transmis¬ 
sion spectroscopy at 461 nm on an atomic beam extracted 
from an oven similar to the one described in section 12.11 
The error signal is obtained using a double pass AOM 
(AOMl in Fig. 12), which is frequency modulated at a 
rate of 10 kHz, with a modulation amplitude of 15MHz. 
The large frequency scan generates an unwanted intensity 
modulation of the spectroscopy beam which is removed by 
a feedback loop on the laser intensity. The transmission 
signal of the laser beam through the atomic beam is col¬ 
lected on a photodetector before being demodulated and 
sent to a PID controller. The feedback loop of the lock 
is closed sending the error signal on the PZT of the mas¬ 
ter laser. We note that the atomic beam has an optical 
thickness of 0.7. The width of the transition is 100 MHz, 
which is Doppler broadened but remains reasonably nar¬ 
row to provide an accurate frequency stabilization. The 
laser beam is perpendicular to the atomic beam to remove 
any bias in the frequency due to Doppler shift. Finally, a 
pair of Helmholtz coils is placed across the atomic beam 
and aligned along the laser beam propagation axis [see 
Fig.[T^a)]. Using a a~ polarized beam, the resonance fre¬ 
quency of the ^Sg,TO = 0—>■ ^Pi,m = —1 transition is 


shifted by Zeeman effect (the Lande factor of the excited 
state is 1). We calibrate this shift as a function of the 
electric current passing through the coils and find the cal¬ 
ibration coefficient to be 36 MHz/A. This scheme allows 
us to scan the lock point of the laser across the strontium 
isotopes as it is illustrated in Fig. [T^b). 



Fig. 12. Optical setup of the 461 nm laser system. The fre¬ 
quencies of the different beams, given on drawing, are relative 
to the lock point. 
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Fig. 13. (a) Schematic representation of the absorption spec¬ 
troscopy system used to lock the frequency of the 461 nm laser 
system. A pair of Helmholtz coils generates a bias magnetic 
field that shifts the atomic transition, (b) The MOT fluores¬ 
cence spectrum for three Sr isotopes. The origin of the x-axis 
corresponds to the ®®Sr blue resonance. The frequency scan is 
performed by changing the value of the bias field. 


3.2 Red laser 

The optical system of the red laser is sketched in Fig. 

We use a master laser made up of an ECDL of 10 cm 
in length. The output coupler is a diffraction grating with 
1800 lines/mm. The facet of the intra-cavity 10 mW diode 
laser is anti-reflection coated. Two isolators with an iso¬ 
lation of 60 dB and 30 dB are placed at the output of 
the ECDL to avoid any optical feedback. The laser is 
frequency locked on the reflection signal of an ultra low 
expansion (ULE) and high finesse {F = 3000) Fabry- 
Perot (FP) cavity using the Pound-Drever-Hall technique. 
The frequency of the master laser is modulated using an 
electro-optic modulator (EOM) at 20 MHz. We generate 
the error signal by demodulating at the same frequency. 
The error signal is sent to the PID controllers which gen¬ 
erate feedback to the current controller and the PZT of 
the master laser for fast and slow corrections respectively. 
The unitary gain of the feedback loop is at 1 MHz. Un¬ 
der this condition, the frequency linewidth is below 1 kHz 
over one second with respect to the FP cavity [H]. The 
FP cavity is isolated from environmental perturbations 
by a vacuum chamber. The vacuum chamber is pumped 
below 10“® mbar, and attached to a double-layer temper¬ 
ature stabilized chamber at 30°C. The long-term drift of 
the FP cavity is below 10 Hz/s. 

The drift of the FP cavity is compensated using satu¬ 
rated fluorescence spectroscopy on the ^Sq —F ^Pi transi¬ 
tion of ®®Sr. The resonance frequency of the ®®Sr is around 
160 MHz below the FP cavity resonance. We use an AOM, 
referred to as AOMl in Fig. in double pass configura- 



®^Sr MOT 
F=9/2^11/2 


'Sr MOT/ 

F=9/2^9/2 


Fig. 14. Optical setup of the 689 nm laser system. 


tion, and driven at an RF frequency of 80 MHz to bring 
the laser beam transmitted through the FP cavity into 
resonance with the atomic transition. We apply a bias 
magnetic field of 10 G to split the magnetic sublevels (m 
= 0, ±1) beyond the Doppler broadening of 9.5 MHz that 
arises from the residual transverse velocity of the atomic 
beam. We perform saturated fluorescence spectroscopy 
on the Zeeman shift free ^So,m = 0 —F ^Pi,to = 0 
transition [Fig[T5|(a)]. The linewidth of the saturated flu¬ 
orescence spectrum is 140 kHz which is limited by the 
mean transit time of the atoms in the beam and inten¬ 
sity broadening [Fig[T5|(b)]. We use a 150 pW probe with 
a beam waist of 10 mm. The probe is retro-reflected us¬ 
ing a cat’s eye setup to minimize the beam misalignment 
which can induce broadening and shifting of the satu¬ 
rated spectroscopy dip. We collect the fluorescence sig¬ 
nal of the atoms on a femtowatt photodetector (Thorlabs 
PDFIOA/M). To increase the fluorescence signal, we also 
collimate and retro-reflect the fluorescence emission in the 


opposite direction, as shown in Fig. 15 a). The signal is 


digitalized and sent to a PC for analysis. The error signal 
is generated by taking the difference of the photodetec¬ 
tor signal while the frequency of the AOMl is changed 
back and forth with a 140 kHz frequency step. We repeat 
this measurement periodically at a rate of 1 s. The error 
signal is sent to a numerical PID controller before be¬ 
ing interpreted in terms of frequency offset. Then, it is 
sent to a Direct Digital Synthesizer (DDS) which controls 
the RF frequency of AOMl. The frequency of the AOMl 
is then broadcasted through a local-area-network system 
and used, in particular, by another DDS to control the RF 
frequency of AOM2. The latter is located at the output 
of the master laser. The accuracy of the frequency sta¬ 
bilization chain is tested on the cold cloud of atoms by 
absorption spectroscopy on the ^So,to = 0 —F ^Pi,to = 0 
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transition. We find a residual frequency bias of 5 kHz and 
a daily variation of ±5 kHz. 

The master laser beam seeds a slave laser, named slavel, 
after being double passed through the AOM2. The out¬ 
put of slavel is further divided into three parts. The first 
part seeds slave2, which is used for the ®®Sr MOT. The 
second part of slavel seeds slaveS which addresses the 
F = 9/2 —>■ F = 9/2 transition used for the ®^Sr MOT 
stirring beam. The third part of slavel passes through an 
EOM at a fixed frequency of 1.2 GHz, generating mainly 
the ±1 lateral orders. We use this beam with multiple 
frequency lines to seed slaved. We tune the parameters 
(temperature and current) of slaved to selectively lock on 
to the -1 order of the frequency spectrum [13]. This laser 
beam is used for the F = 9/2—= 11/2 MOT transi¬ 
tion of ®’’Sr. 




Fig. 15. (a) Optical setup used for the saturated fluorescence 
spectroscopy on the ^So —Pi transition of ®®Sr. (b) Typical 
saturated fluorescence spectroscopy signal. The blue curve is 
the experimental spectrum whereas the red curve is a fit of the 
experimental signal summing a broad gaussian and a narrow 
lorentzian shape. 


3.3 Repumper lasers 

Due to the presence of the low lying energy state, the 
blue cooling transition is not completely closed. This leads 
to the shelving of atoms in the ®P 2 state (see Fig. 0.TO 
optically repump the ®P 2 population to the ground state, 
we drive the ®P 2 —®Si transition at 707 nm. Since some 
atoms can be repumped to the long lived ®Po state as 


well, we also drive the ®Po —>■ ®Si transition at 679 nm. 
The 679 nm and 707 nm lasers are home-made ECDLs. 
The frequencies of these lasers are monitored and locked 
by digital feedback on a wavelength meter (HighFinesse, 
WS/7). The accuracy in frequency of the wavelength me¬ 
ter is dO MHz. It is sufficient since the efficiency of re¬ 
pumping seems to be not affected over a frequency range 
of ~ 100 MHz across the resonance, most likely because 
of the spatial inhomogeneities in the transition frequency 
induced by the MOT magnetic field gradient. 



Fig. 16. Repumping scheme for ®®Sr and ®^Sr. (a) The black 
solid lines are the different transition for ®^Sr and ®®Sr. The 
blue curves are the laser spectra after passing through the 
EOM, modulated at 540 MHz. Those spectra are revealed us¬ 
ing a low finesse {F = 100) FP frequency analyzer. The num¬ 
bers, ranging from —3 to -|-3 at the vicinity of the peaks, in¬ 
dicate the lateral bands order of the EOM. (b) Fluorescence 
signal of the ®®Sr and ®^Sr MOT as a function of the 707 nm 
repumper optical power. The 679 nm laser is maintained at its 
maximum power, (c) Fluorescence signal of the ®®Sr and ®^Sr 
MOT as a function of the 679 nm repumper optical power. The 
707 nm laser is maintained at its maximum power. 


The bosonic case is rather simple, since only two tran¬ 
sitions have to be considered. For ®^Sr, the situation is 
more complicated because of the large hyperfine structure 
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of the ^P 2 and ^Si levels. The transitions ^P 2 —i ^Si and 
^Pg —> ^Si are shown in Fig. 16 a). In order to pump 
the population in the ^Pg level, we need only one transi¬ 
tion as there are no hyperfine levels in this state. 679 nm 
repumper laser addresses the F = 9/2—i-F = 9/2 transi¬ 
tion, which is 550 MHz away from the ®®Sr transition. An 
efficient repumping scheme at 707 nm is more challeng¬ 
ing, as we need to repump five hyperfine levels, namely 
F = 5/2, 7/2,9/2,11/2, and 13/2. We choose to address 
the transitions F = 5/2 F = 7j2, F = 7j2 ^ F = 9/2, 
F = 9/2 -)■ F = 9/2, F = 11/2 F = 9/2, and F = 
13/2 —>• F = 11/2, because these transitions are closer 
to ®®Sr repumping line and they are almost equally sepa¬ 
rated. We generate a multiline laser using an EOM driven 
at 540 MHz for the 707 nm and 679 nm laser beams. Both 
laser beams are mixed using a polarizing beam splitter, 
thus they have orthogonal polarization. The polarization 
of the 707 nm laser is chosen such that the efficiency of 
the EOM is optimum. The modulation index of the EOM 
is adjusted to 1.8 to get substantial power into the central 
band and in the ±1 and ±2 lateral bands [see Fig. [T^a)]. 
The efficiency of the EOM for the 679 nm laser beam is 
much weaker but it is sufficient to repump the atoms from 
the ^Pg state. The fluorescence signal from the cold cloud 
in the blue MOT is shown in Fig. 16 [b,c) as a function of 
the 679 nm and the 707 nm laser power. With this scheme, 
we are able to repump ®®Sr and ®^Sr without changing the 
frequency lock points of the lasers. 


4 Conclusion 

We described a new experimental apparatus that produces 
ultracold gases of ®®Sr and ®’^Sr. The system is character¬ 
ized by a high loading rate of the blue MOT done on the 
dipole allowed transition. This loading is performed under 
high vacuum environment. We achieve this by introduc¬ 
ing an atomic beam deflector (2D MOT) after the Zee- 
man slower. Using the 2D MOT, we selectively deflect the 
longitudinally decelerated atomic beam into the science 
chamber. We show that the 2D MOT efficiency is limited 
by the decay into the ^D 2 long-lived state. The maximum 
loading rate measured in the MOT is 6 x 10® s“^. A dif¬ 
ferential pumping tube is place between the 2D MOT and 
the science chamber to reduce the background pressure at 
the MOT location. The lifetime of the atoms in the mag¬ 
netic trap is around 54 s. An extra optical axis is opened 
up along the atomic beam direction after this deflection. 
This can be useful, for example, to further manipulate the 
beam towards a continuous atom laser. 

The frequency lock point of the 461 nm laser system 
can be tuned by shifting the atomic transition with Zee- 
man effect. This technique allows us to easily tune the 
laser on the desired isotope. The loss mechanism induced 
by spontaneous decay of the excited state of the cooling 
transition ^Pi to the underlying long live state ^F >2 are 
properly shelved using a simple multi-line 707 nm laser. 
Finally, we show that a transfer efficiency of 50% from 
the blue MOT to the narrow intercombination line red 
MOT is achieved. 


5 Acknowledgement 

This work was supported by the CQT/MoE funding Grant 
No. R-710-002-016-271. 


References 

1. H. Katori, T. Ido, Y. Isoya, M. Kuwata-Gonokami, Phys. 
Rev. Lett. 82, 1116 (1999) 

2. T. Mukaiyama, H. Katori, T. Ido, Y. Li, M. Kuwata- 
Gonokami, Phys. Rev. Lett. 90(11), 113002 (2003) 

3. F. Sorrentino, G. Ferrari, N. Poli, R. Drullinger, G.M. 
Tino, Modern Physics Letters B 20(21), 1287 (2006) 

4. W. Shao-Kai, W. Qiang, L. Yi-Ge, W. Min-Ming, L. Bai- 
Ke, Z. Er-Jun, L. Tian-Chu, F. Zhan-Jun, Chinese Physics 
Letters 26(9), 093202 (2009) 

5. S. Stellmer, M.K. Tey, B. Huang, R. Grimm, F. Schreck, 
Phys. Rev. Lett. 103, 200401 (2009) 

6. Y.N.M. de Escobar, P.G. Mickelson, M. Yan, B.J. DeSalvo, 

S. B. Nagel, T.G. Killian, Phys. Rev. Lett. 103, 200402 
(2009) 

7. M. Chalony, A. Kastberg, B. Klappauf, D. Wilkowski, 
Phys. Rev. Lett. 107, 243002 (2011) 

8. S. Stellmer, B. Pasquiou, R. Grimm, E. Schreck, Physical 
review letters 110(26), 263003 (2013) 

9. K.Y. Khabarova, A.A. Galyshev, S.A. Strelkin, A.S. 
Kostin, G.S. Belotelov, O.I. Berdasov, A. Gribov, N.N. 
Kolachevsky, S.N. Slyusarev, Quantum Electronics 45(2), 
166 (2015) 

10. M. Takamoto, H. Katori, Phys. Rev. Lett. 91(22), 223001 
(2003) 

11. M.M. Boyd, A.D. Ludlow, S. Blatt, S.M. Eoreman, T. Ido, 

T. Zelevinsky, J. Ye, Phys. Rev. Lett. 98(8), 083002 (2007) 

12. A. Derevianko, H. Katori, Rev. Mod. Phys. 83, 331 (2011) 

13. A.D. Ludlow, M.M. Boyd, J. Ye, E. Peik, P.O. Schmidt, 
Rev. Mod. Phys. 87, 637 (2015) 

14. N. Poli, C.W. Oates, P. Gill, G.M. Tino, Rivista Nuovo 
Cimento 12, 555 (2013) 

15. B. Bloom, T. Nicholson, J. Williams, S. Gampbell, 
M. Bishof, X. Zhang, W. Zhang, S. Bromley, J. Ye, Na¬ 
ture 506(7486), 71 (2014) 

16. 1. Ushijima, M. Takamoto, M. Das, T. Ohkubo, H. Katori, 
Nature Photonics 9(3), 185 (2015) 

17. C.W. Chou, D.B. Hume, T. Rosenband, D.J. Wineland, 
Science 329(5999), 1630 (2010) 

18. T. Rosenband, D.B. Hume, P.O. Schmidt, C.W. Chou, 
A. Brusch, L. Lorini, W.H. Oskay, R.E. Drullinger, T.M. 
Fortier, J.E. Stalnaker et ah. Science 319(5871), 1808 
(2008) 

19. D. Akamatsu, M. Yasuda, H. Inaba, K. Hosaka, T. Tanabe, 
A. Onae, F.L. Hong, Optics express 22(7), 7898 (2014) 

20. G. Tarallo, M. T. Mazzoni, N. Poli, V. Sutyrin, D. 
X. Zhang, M. Tino, G. Phys. Rev. Lett. 113, 023005 
(2014) 

21. B.J. DeSalvo, M. Yan, P.G. Mickelson, Y.N. Martinez de 
Escobar, T.C. Killian, Phys. Rev. Lett. 105(3), 030402 
( 2010 ) 

22. P.G. Mickelson, Y.N. Martinez de Escobar, M. Yan, B.J. 
DeSalvo, T.G. Killian, Phys. Rev. A 81, 051601 (2010) 

23. S. Stellmer, M.K. Tey, R. Grimm, F. Schreck, Phys. Rev. 
A 82, 041602 (2010) 




12 


T. Yang et al.: A high flux source of cold strontium atoms 


24. M.K. Tey, S. Stellmer, R. Grimm, F. Schreck, Phys. Rev. 
A 82, 011608 (2010) 

25. M. Foss-Feig, M. Hermele, A.M. Rey, Phys. Rev. A 81, 
051603 (2010) 

26. X. Zhang, M. Bishof, S. Bromley, C. Kraus, M. Safronova, 
P. Zoller, A. Rey, J. Ye, science 345(6203), 1467 (2014) 

27. D. Banerjee, M. Bogli, M. Dalmonte, E. Rico, P. Stebler, 
U.J. Wiese, P. Zoller, Phys. Rev. Lett. 110, 125303 (2013) 

28. Y.X. Hu, C. Miniatura, D. Wilkowski, B. Gremaud, Phys. 
Rev. A 90, 023601 (2014) 

29. Y. Bidel, B. Klappauf, J.G. Bernard, D. Delande, 
G. Labeyrie, G. Miniatura, D. Wilkowski, R. Kaiser, Phys. 
Rev. Lett. 88, 203902 (2002) 

30. T. Chaneliere, D. Wilkowski, Y. Bidel, R. Kaiser, 
G. Miniatura, Phys. Rev. E 70, 036602 (2004) 

31. M. Chalony, R. Pierrat, D. Delande, D. Wilkowski, Phys. 
Rev. A 84, 011401 (2011) 

32. G.C. Kwong, T. Yang, M.S. Pramod, K. Pandey, D. De¬ 
lande, R. Pierrat, D. Wilkowski, Phys. Rev. Lett. 113, 
223601 (2014) 

33. J.M. Danet, M. Lours, E. De Glercq et ah, Ultrasonics, 
Ferroelectrics and Frequency Control, IEEE Transactions 
on 61(4), 567 (2014) 

34. S. Stellmer, R. Grimm, F. Schreck, Phys. Rev. A 87(1), 
013611 (2013) 

35. G. Labeyrie, E. Vaujour, C.A. Mueller, D. Delande, 

C. Miniatura, D. Wilkowski, R. Kaiser, Physical review 
letters 91(22), 223904 (2003) 

36. H. Pauly, Atom, Molecule, and Cluster Beams I: Cluster 
Beams, Fast and Slow Beams, Accessory Equipment and 
Applications, Vol. 1 (Springer Science & Business Media, 
2000) 

37. P. Berthoud, A. Joyet, G. Dudle, N. Sagna, P. Thomann, 
EPL 41(2), 141 (1998), ISSN 0295-5075 

38. H. Chen, E. Riis, Applied Physics B 70(5), 665 (2000) 

39. K. Dieckmann, R.J.C. Spreeuw, M. Weidemiiller, J.T.M. 
Walraven, Physical Review A 58(5), 3891 (1998) 

40. T.P. Dinneen, K.R. Vogel, E. Arimondo, J.L. Hall, A. Gal¬ 
lagher, Physical Review A 59(2), 1216 (1999) 

41. A. Caires, G. Telles, M. Mancini, L. Marcassa, V. Bag- 
nato, D. Wilkowski, R. Kaiser, Brazilian journal of physics 
34(4A), 1504 (2004) 

42. M. Yasuda, H. Katori, Phys. Rev. Lett. 92, 153004 (2004) 

43. S. Stellmer, F. Schreck, Physical Review A 90(2), 022512 
(2014) 

44. N. Poll, R.E. Drullinger, G. Ferrari, J. Leonard, F. Sor- 
rentino, G.M. Tino, Phys. Rev. A 71, 061403(R) (2005) 

45. K.R. Vogel, T.P. Dinneen, A. Gallagher, J.L. Hall, IEEE 
Transactions on Instrumentation and Measurement 48, 
618 (1999) 

46. T. Chaneliere, J.L. Meunier, R. Kaiser, C. Miniatura, 

D. Wilkowski, Journal Of The Optical Society Of America 
B Optical Physics 22(9), 1819 (2005) 

47. T. Chaneliere, L. He, R. Kaiser, D. Wilkowski, The Euro¬ 
pean Physical Journal D 46, 507 (2007), ISSN 1434-6060, 
1434-6079 

48. Y. Bidel, PhD thesis, unpublished (2004) 

49. M.S. Pramod, T. Yang, K. Pandey, M. Giudici, 
D. Wilkowski, Eur. Phys. J. D 68(7), 1 (2014), ISSN 1434- 
6060, 1434-6079 



